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Abstract — Of concern in this work is the derivation and implementation ofthe multistep methods through Taylor’s expansion 
and numerical integration. For the Taylor’s expansion method, the series is truncated after some terms to give the needed 
approximations which allows for the necessary substitutions for the derivatives to be evaluated on the differential equations. 
For the numerical integration technique, an interpolating polynomial that is determined by some data points replaces the 
differential equation function and it is integrated over a specified interval. The methods show that they are only convergent if 
and only if they are consistent and stable. In our numerical examples, the methods are applied on non-stiff initial value 
problems of first-order ordinary differential equations, where it is established that the multistep methods show superiority 
over the single-step methods in terms of robustness, efficiency, stability and accuracy, the only setback being that the multi- 
step methods require more computational effort than the single-step methods. 
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I. INTRODUCTION 

Linear multistep methods (LMMs) are very popular for solving initial value problems (TVPs) of ordinary differential 
equations (ODEs). They are also applied to solve higher order ODEs. LMMs are not self-starting hence, need starting values 
from single-step methods like Euler’s method and Runge-Kutta family of methods. 

The general fc-step LMM is as given by Lambert [1] 

k k 

'Z l ct j x„ + j = hY J Pj<p n+ j ( 1 ) 

j=0 j=0 

where a . and ff are uniquely determined and a 0 + ff, ^ 0, a k = / . The LMM in Equation (1) generates discrete schemes 

which are used to solve first-order ODEs. Other researchers have introduced the continuous LMM using the continuous 
collocation and interpolation approach leading to the development of the continuous LMMs of the form 

y(0=Yj a j( 1 K +j = 1 ton+s (2) 

]=0 j=0 

where 0 Cj and ff j are expressed as continuous functions of x and are at least differentiable once [2]. 

According to [3], the existing methods of deriving the LMMs in discrete form include the interpolation approach, numerical 
integration, Taylor series expansion and through the determination of the order of LMM. Continuous collocation and 
interpolation technique are also used forthe derivation of LMMs, block methods and hybrid methods. 

In this study, we present the general multistep method, some of its different types and examine their characteristics, hi light 
of this, we investigate the stability and convergence of these methods, compare the multistep methods with the single-step 
methods in operational time, accuracy and user-friendliness via some numerical examples. 

hi practice, only a few of the initial value differential equations that originate fromthe study of physical phenomena have 
exact solutions. The introduction however, of the multistep methods as numerical techniques is used in finding solutions to 
problems that have known exact solutions and in extension handle those problems whose exact solutions are not known. We 
shall limit this study to only non-stiff initial value problems of first-order ordinary differential equations. 
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The Linear Multistep Methods 

The general linear multistep method is given by 

k k 

Y,cc j x n+j =hY j P j (p n+j (3) 

j=0 j=0 

where the and /T are constants, whereas a k =0 and both a 0 and fi 0 are not zero. Since (3) can be multiplied on both 
sides by the same constant without altering the relationship, the coefficients a. and /?. are taken arbitrarily to the extent of a 
constant multiplier. In this work however, we will assume that a k = 1 ■ If J3 k =0 , the equation (3) is explicit otherwise it is 
implicit [4]. 

The Adams methods 

These are the most important linear multistep methods for non-stiff initial value problems. It is the class of multistep 
methods (3) with a k = 1 , a k _ t = -1 and = 0 , j =0,1,2, ■■■,k —2 . If equation (3) is given by 

k-l 

X n+k = X a+k _, +hY,Pj(p n+j , (4) 

i=o 

then we have the Adams-Bashforth methods. And, if it is given by 

k 

X n+k =x n+k _ 1 +h£ d 0 j <p n+j , (5) 

i=o 

then we have the Adams-Moulton methods [5]. 

Predictor-Corrector (P-C) method 

The multistep methods are often implemented in a ‘predictor-corrector’ fomi In this way, a preliminary calculation is done 
using the explicit form of the multistep method then corrected using the implicit form of the multistep method. This is done 
by two calculations of the function rp at each step of this computation. 

Order of linear multistep methods 

We can associate the linear multistep method (3) with the linear difference operator 0 , defined by 

0[x(r );h\ = Yj[ a jX( r + jh )-hPjX'( r + jhf\ , (6) 

j=0 

where x ( r ) is any arbitrary function that is continuously differentiable on the interval [a,b]. If the operator is allowed to 
operate on an arbitrary test function x(r) with as many higher derivatives as we require, we can formally define the order of 
the operator and of the associated multistep method without invoking the solution of the initial value problem (4) which may 
possess only a first derivative. If we expand the test function x( r + jh ) and its derivative x'( r + jh) as Taylor series about r 
and collect terms in (6) we have 

0[x( r );h] = C 0 x( r ) + Cjhx'( r)- 1-f C q h q x (q) ( r)-\ — (7) 

where the C q constants [1]. 


Definition 1 The difference operator (6) and the linear multistep method (3) associated with it are of order 8 if, in (7), 
C 0 = C t =C 2 = ■■■ = C s = 0 and C s+I ^0 . 

The following formulae for the constants C q in terms of the coefficients (Xj and /T are given as: 

C 0 = a 0 + a, + a 2 H- 1 -a k (8) 

Cj = cCj + 2a 2 + 3a 3 -4-f ka k —(/3 0 + /3 k + /3 2 -\ - P k ) (9) 


C = — (a, +2 q a , H- 1 - k q a k - 


(q—i)!(Pi 


-2 q -‘ p 2 +--- + k q ~ 1 P k ) 


for q =2,3,-■ ■ [1]. 


( 10 ) 
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II. CHARACTERISTICS OF THE METHODS 

With the number of approximations involved during computations using the multistep methods, the problem of consistency, 
stability and convergence call for discussion. The approximation in a one-step method depends directly on previous 
approximations alone, while the multistep method uses at least two of the previous approximations. 

Consistency 

The linear multistep method (3) is consistent if it has order 5>1[ 4]. From (8), (9) and (10), it follows that the method is 
consistent if and only if the following two conditions hold. 


k 


Z a j= 0 ’ 

i=o 

(11) 

k k 

Z;'« 7 =Z^ 

j=0 j=0 

(12) 


[1]. 

We shall subsequently considerthe limits as /r—>0, n—> co and nh—>r-a remaining fixed. 

Let x n tend to x( r ) in the limit, that is x n —>x(r). Since k is fixed, we have that x n+j —> x( r ), for j =0,1,2,-■■ ,k , or 
x(r) = x n+j + Tj J h), j =0,1,2,---,k 
where lim i v Jh) = 0 , j = 0,1,2, ---,k [1]. 

Hence, we have 

k k k 

Yj a j X ( r >= H a j X n+i +X a j T jJ h > ’ 

j=0 j=0 j=0 

or replacing the first term on the right hand side of the equation by the temion the right hand side of (3) we have 

k k k 

x( rj^cij = + 'Z a Pj,( h). 

j=0 j=0 j=0 

In the limit, both temis on the right hand side vanish. Therefore, the left hand side becomes zero. The left hand side is not in 
general equal to zero, so we conclude that y a. = 0. The above argument holds if we merely assume that {-*(/)} tends to 

some function x(r). 

Condition (12) ensures that the function x(r) does in fact satisfy the differential equation. For, under the limiting process, 

——- -—>x'(r), for j = l,2,---,k, 

jh 

or, x n+j - x n = jhx'( r) + jh<p Jn (h), for j = l,2,---,k, 
where lim (p j n { h ) = 0 . Hence, 

k k k k 

Y, a i X n + j -Tk a i X n = h Tj-> a j X '( r ) + H a i T jJ h > 

j=0 j=0 j=0 j=0 

or 

k k k k 

bYsPjVn+j - X nH, a j =hX '( r >TjJ a i +h Yjj <X i < PiJ k ) • 

j=0 j-0 j-0 j=0 

Since ^ n0 a j - 0 , we have, on dividing through by h , 

k k k 

X Pj<Pn + j = X '( r + Z j a i ( Pj.n ( h >-- 

i=o j=o j=o 

Under the limiting process, <p n+j —> <p( r,x( r )), and, in the limit, 

k k 

<P(r,x(r ))J] Pj = x'( r ja j .. 

i=o j=o 

k k 

Thus x( r) satisfies the differential equation (4) if and only if ^ - 0 i a j = Z Pi ■ shows that if the sequence 

j=o 

converges to the solution of the initial value problem (4) then the conditions (11) and (12) must hold [1]. 
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Stability 

Let us introduce the first and second characteristic polynomials of the multistep method (3), defined as p( A ) and <j(A) 

k k 

respectively, where p( A ) = 2_ j a j A J and a( A ) = y' j8 i A i 

j=o j=0 

[ 6 ]. 

It follows from conditions (11) and (12) that the linear multistep method is consistent if and only if p(l) = 0 and 
p'( 1 ) = <j( 1) [1]. The stability of the multistep technique with respect to round-off error is clearly dictated by the magnitude 
of the zeros of the first polynomial above. However, the methods we have discussed in this work are zero -stable by virtue of 
their characteristics. The following are motivated by the types of zeros of the characteristic polynomial. 

Theorem 3.1 

Let Aj , A 2 , ..., A k denote the roots (real or complex), which may not necessarily be distinct, of the characteristic equation 
associated with the multistep difference method. If \A n \ < 1 , for n = l,2,---,k and all the roots with absolute value equal to 1 
are simple roots, then the difference method satisfies the root condition [7]. 

Theorem 3.2 

i) The methods that satisfy the root condition with A = 1 as the only root of the characteristic equation with magnitude 
equal to 1 are said to be strongly stable. That is, the roots lie on the unit disc. 

ii) If a method satisfies the root condition and has more than one distinct root with magnitude equal to 1, it is said to be 
weakly stable. 

iii) If a method does not satisfy the root condition, it is said to be unstable. A multistep method is said to be stable if and 
only if it satisfies the root condition [7], 

Convergence 

One basic property that is demanded of an acceptable linear multistep technique is the convergence of the solution that 
is generated by the method, in some sense, to the theoretical solution x(rj as the step-size h goes to zero. A linear multistep 

method is convergent if and only if it is consistent and stable, otherwise it is not convergent [8]. If a method is consisten t but 
not stable, then it is not convergent. Also, if a method is stable but not consistent then it is not convergent. 

Obtaining Starting Values 

A multistep method is not self-starting, that is, a k-step multistep scheme requires some k previous values x 0 ,x I ,x 2 ,...x k _ I . . 
These k values that are needed to start the application of the multistep method are gotten by a single step method such as 
Taylor series method, Euler method or Runge-Kutta method. The starting method should be of the same or even lower order 
than the order of the multistep method itself. 


Taylor series method 

Let us consider the initial value problem 

x' = cp(r,x), x 0 =a. (13) 

Let us consider a numerical solution to (13) above using a k-step multistep method of order 5 . We require that the starting 
values x ( ., i = 1,2,- ■ -,k —1 should be calculated to an accuracy that is at least as high as the accuracy of the multistep method 


itself. That is, we require that x t — x( r i ) = 0(li s+I ), i = 1,2,-- -,k —1 . 

If enough partial derivatives of cp( r,x) with respect to r and x exist, then we will use a tmneated Taylor series to estimate 
x i to any required degree of accuracy [9]. Thus, we have 


1 2 h 3 

x !+1 = x( r t )+ hx'( r t ) + x"( r { )+—x m ( r t ) -4- bO( h S+I ) , 


(14) 


for i = l,2,---,k—l. The derivatives in (14) are evaluated by successively differentiating the differential equation. Thus, 
x(r 0 ) = x 0 , 
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x'( r) = cp{ r,x), 


, dtp dtp , drp d(p 

dr 8x dr dx 

„ , d 2 m „ d 2 cp , d 2 m drpdrp 

dr' drdx dx drdx 


d<P V 
dx J 


This approach is theoretically flawless. Nevertheless, the evaluations of the total derivatives can be excessively tedious and 
may not be adopted for an efficient computation. 


Euler method 

This is another method that can be employed to generate all the needed starting values for a linear multistep method. 

Consider the equation 

x’ - rp( r, x ), a < r < b , x( a ) = a. 

Let us suppose that the solution to the initial value problem (15) above has two continuous derivatives on the interval , 
so that for each i = 7,2, •••,1V —1 , we have 

x M =x(r i ) + ( r M - r. )x'( r { ) + ^ r ' * x’( £. ), 

for some g. in (r,r + ,) .Since h = r i+/ — r , we have 

x i+I =x(r) + hx'( r i) + Y, x "( & ) 

and, since x(r) satisfies our differential equation, we have 

h 2 

x i+I = x( r { ) + h<p( r,x( r t )+ — x "(^ ) . 

By deleting the remainder term, the Euler method becomes 
x 0 = a 

x j+ i = x { + h<p( r, Xj ), 
for each i = 1,2,-■■ ,N -1. 

The Euler method is gotten when the Taylor series method above is of order 8 = 1. The simplicity of this method may be 
used to illustrate the techniques we intend to adopt in starting the multistep methods. 

Runge-Kutta method 

This method can also be applied to generate starting values for any multistep method. 

We consider the equation 

M 

x( r + h ) = x( r ) + h^ wk : (15) 

S=1 

where 

s—l 

k I =cp(r,x), k s =<p( r + a s h,x+h^ j j3 si k i , for s = [1], 

1=7 

We call this an M -order Runge-Kutta method and it involves M function evaluations at each step. Each k, , k =1,2,---,M , 
may be interpreted as an approximation to the derivative x!(r). 

The objective is to choose w s , a s and J3 sr so that the coefficients of H , i =1,2,---,M , in equation (14) are identical with 
those of the equation (15). That is, the method must compare with the Taylor series method after its expansion. The higher 
order derivatives of the Taylor series expansion is given by 
x' = <p 

x" = <p r + <p x x' = (p r + <p<p x 

X m = <Prr + 2( P<P rI + <P 2 <P a + (P x (Pr + W’ 
and so on. 
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We note that the Runge-Kutta methods are not unique due to the manner in which they are derived. However, any Runge- 
Kutta methods of the same order are equivalent. 


Runge-Kutta method of order two 

This method uses two evaluations and it is given by 
x i+1 = Xj + hw.k i + hw 2 k 2 , 

where k, = cp( r t ,x j ), k 2 =rp( r i +^h,x i +^hk 1 ) 

Runge-Kutta method of order four 

This method uses four evaluations. It is given by 

Xj + j — Xj H— (k, + 2k 2 + 2k 2 + k 4 , 

6 


( 16 ) 


Where k, = h(p{ r i ,x j ) 

k , = ha/ r. + — h,x. H— hk,) 

2 , 2 '2 

h h 

k 3 =htp(r i +-,x i +-k 2 ) 
k 4 = hcp( +h, x i + k 3 ) 

[ 10 ]. 

We note however, that this is not unique. The Runge-Kutta method of order four shall be used in obtaining the starting values 
for the implementation of the multistep methods adopted in this work. 


Derivations 

Any specific linear multistep may be derived in a number of different ways. We shall consider a selection of different 
approaches which cast some light on the nature of the approximation involved. 


Derivation through Taylor expansions 
Euler method 

Let us consider the Taylor series expansion for x( r n +h) about r n , 


x( r n + h) = x( r n ) + hx'( r ) + 


/r 

2! 


x’( r n )-\— 


(17) 


If we tmncate this expansion after two terms and substitute for x'( r) from the differential equation (4), we have 
x( r n + h ) = x( r n ) + hcp( r n , x( r n )) , 
h 2 „ „ 

and the tmncation error is — x ( q n ) for <^ n in (r n ,r n+1 ). 

If x(r n ) and x( r +h) are replaced by x n and x n+1 , we get 

X n + l= X „+ h <P„’ ( 18 ) 

which is an explicit linear one-step method [11]. This shall be used in solving the numerical examples in this work. 


Mid-point rule 

Let us considerthe Taylor series expansions for x( r n +li) and x( r n -h) about r n Thus, 

f h 2 „ h 3 

x( r n + h ) = x( r n ) + hx (r n )+—x (r n ) +—x ( r n ) H— 

,2 ,3 

x( r n - h ) = x( r n ) - hx'( r n ^ + ~^j x "( r » ) _ “7 xm ( r „ )3 — 

Subtracting, we have 

h 3 

x( r n + h) — x( r n - h) = 2hx'( r n ) + — x"'( r n ) H— 

Replacing x( r n + h ) and x( r n -h) by x n+1 and x n l , we have x n+I -x n _, = 2hcp n . 
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This can be brought into the standard formof the linear multistep method (13), after replacing n by n + 1, as 

x n + 2~ x n = 2h( P„ + , (! 9 ) 

Its tmncation error is ±—h 3 x m (£ ) [1]. 

3 

This is the mid-point mle and it shall be used in solving the numerical examples presented in this work. 


The trapezoidal rule 

If we wish to find the most accurate one-step implicit method x n+I +a 0 x n = h( P I (p n+I + P,/P n ), we write down the associated 
approximate relationship 

x (r n +h) + a 0 x(r n )*h[P,x'( r n +h) + P 0 x '(r j] (20) 

and choose a 0 , P 0 and p so as to make the approximation as accurate as possible. 

The following expansions are used: 

, h 2 „ 

x( r + h ) = x( r ) + hx ( r ) H- x ( r ) H— 

" 2 / 

h 2 

x'( r + h ) = x'( r ) + hx"( r ) H- x"'( r ) -t— 

" 2! n 

Substituting these two equations into (20) and collecting the terms on the left-hand side gives 
C 0 x( r n ) + C,lix'( r n ) + C 2 h 2 x”( r n ) + C 3 h 3 x”'( r n ) H— « 0 , 

where 


Cg — 1 + Cl g, 

c, = i-p,-p 0 , 

c -r A - 


Thus, in order to make the approximation in equation (20) as accurate as possible, we choose a 0 =—l, P s = P 0 = —. The 

value of Cj then becomes ~~ ■ The linear multistep method is now 

h 

X n + 1 +X n =^(<Pn + l+<Pn)- 

This is the trapezoidal rule and its local tmncation error is ± — h 3 x'”( E, n ) [1]. 


Derivation through numerical integration 

This technique can be used to derive only a subclass of linear multistep methods consisting of those methods for which 
a k = +1 , cCj = —1 , or. = 0 , i =0,1,2, ■■■,_/ —l,j +1 , j ;/= k . To start the derivation of any multistep method, we should note 

that the solution of the initial value problem given as 
x' = cp{ r,x ), a<r <b , x( a ) = a , 
if integrated over [ r n , r n+I ], has the property that 

x( r n+I )-x( r n ) = J " +( x'( r )dr = J cp( r,x( r ))dr . 

Consequently, 

x( r n+1 )-x(r n ) = J '”"' rp(r,x(r ))dr . (21) 

We will integrate some interpolating polynomial P(r) to <p( r,x( r )) which is determined by some previous data points that 
were obtained. Then, equation (21) becomes 

x( r n+I )-x(r n ) = J r " + ' P( r )dr . (22) 
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Simpson’s rule 

Suppose we want to derive a two-step method, we consider the identity 

x( r n 2 ) — x(r n ) = J" + “ x'( r )dr . (23) 

Using the differential equation (4), we can replace x by tp( r,x) in the integrand. The only available data for the 
approximate evaluation of the integral will be the values cp n , rp n+1 and rp n+2 . Let P(r) be the unique polynomial of second 
degree passing through the three points ( r n ,rp n ), (r n+I ,rp +/ ) and ( r rl+2 ,(p rl _,). By the Newton-Gregory forward interpolating 
formula, 

j 1 / j 1 _ J \ 

P(r) = P( r n +sh) = (p n + sArp n + ——— A 2 <p n . 

We now make the approximation 

f ” 2 x'( r )dr « f P( r )dr = f [<p + sArp +-s(s-l)A 2 (p ]hds 

J r„ J r n JO 2 

= h[2<p n +2Acp n +^A 2 (p n ]. 

If we expand Arp n and A 2 cp n in terms of <p n , (p n+I , (p n+2 and substitute in (23), we have 
h 

* n+ 2 -*„=-[ <Pn + 2 + 4( Pn+l + <P,J ' 

Then, the truncation error becomes ± — h 5 x <5) ( £ ) . 

90 " 

This is the Simpson’s mle and it is the most accurate implicit linear two-step method [1]. 
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Adams-Bashforth methods 

Though any fonn of the interpolating polynomials could be used for the derivations, the Newton backward-difference 
formula will be used for the purpose of convenience. For us to derive an explicit A:-step Adams-Bashforth method, we now 
form the backward-difference polynomial P k _,( r) through (r n ,tp(r n ,x(r n ))), (r n I ,cp( r n I ,x(r n _ I ))) , 

( r n+l _ k ,qir n+l k ,x(r n+l k ))) . Since P k I (r) is an interpolating polynomial of degree k-1 , then for some £ n in ( r n+I k ,r n ), 
we have 

rp <k> ( £ ,x( £ )) 

rp( r,x( r) = P k _ 1 (r) + - - y- -— ( r-r n )( r-r n _ 2 )■■■( r-r ,_ k ) . 

Ic ! 


Introducing the substitution r = r n + sh , and with dr = hds into P k l and with the error term, it implies that 


<P(r,x(r ))dr = | ^( -1)" 


V '"<P( r „’ x ( r „ ))dr 


4 : 


^<p ,kl (£ n ,x(£ n )) 


k! 


(r-r„)(r- r n _ 2 )■ ■ ■ ( r - r n+I _ k )dr 


fy m < P (r n ,x(r n m-l) m \ , \ ^ds+’^-l' S (s + l)-(s + k-l)<p (k, (£ n ,x(£ n ))ds . 
TTn J0 m k J0 


The integrals (-1 ) m 


f' 

M 

Jo 



ds 


for various values of in can be evaluated easily as displayed below. 
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As a consequence of these evaluations, 
j <p( r,x( r ))dr = h 

4— -h k+1 j o s( s + 1 )■■■(s + l-k )rp' k Y <^ n ,x())ds . 


tp( r n ,x( r n ))-\~V(p( r n ,x( r n ))+ ~^V 2 <p( r n ,x( r n )) + ■ 
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(24) 


The coefficient s( s + 1)(s + 2 )■■■( s + k-1) does not change sign on the interval \0,l\ [12]. The weighted mean value 
theorem for integrals can be used to deduce that for a number ju n , where r n+I _. < fi rl < r n+I , the error term in equation (24) 
becomes 


~rr\ n s ( s + 1 )'"( s + 1 ~k)<P ik> (L’X())ds = -— s (s + 1 )■■■(s + k-1 )ds 
k! J0 k! J0 

or the error term can be written as 

h k+1 <p <k> (ju n ,x( J u n ))(-l) m \ I 0 \ : ds 

in 


r' 


Jo 

U J 


Adams-Bashforth two-step explicit method: 

X ( r n + l)~ X ( r n ) + M ^ r n’ X ( ))+^<P( T n> X ( ))1 

h 

= x( r )+hcp{ r n ,x( r )) + -[<# r n ,x( r n ))-rp( r n _,,x( r n _, ))\ 
h 

= x(r„)+-[3<p( r n ,x( r n ))-<p( r n _ t ,x( r n _,))]. 

Consequently, the explicit Adams-Bashforth two-step method is: 
x 0 = a 0 , Xj = a } 

h 

X n + 1 = X n +~[3<P( r nM ^ ))~ <P( r „-l’ X i r „-l ))l ■ 

This can be taken into the standard form of the linear multistep method, after replacing n by n + 1. Thus, 


X n + 2= X n + l +, ^[ 3( Pn + l-(Pn l 
The local truncation error is 

T nJ h ) = J 2 x '( M n )h 2 
for some ju n e (r n _ 17 r n+I )[ 7], 

This method shall be used in solving the numerical examples in this work. 


(25) 


Adams-Bashforth three-step explicit method: 

x(r ii+1 )-x(rj~h 


(Pi r „’ x ( r „ ))+ 1 ~^(Pi r „ ’ x i r „ )) + j 2 ^ 2< Pi r n ,x( r )) 


= x( r n )+h 


(Pi r „’ X ( r n ))+~{(pi r n’ X ( r n ))~(Pi r n-lM Y n-1 ))} 


+ — {(Pi r n’ x i r „ ))-2<P( r n -„x( r n _j ))+cp( r n _ 2 ,x( r n _ 2 ))} 


= x ( r n )+—[ 23( Pi r n’ X i r n ))~^(pi L-n X ( L-l )) + 5( Pi r „-2’ X i V n-2 ))]■ 
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Consequently, the Adams-Bashforth three-step explicit method is 
x 0 = a 0 , x 1 =a r x 2 = a 2 

X n + 1 = X n + Jj[ 23< ti r n’ X n )~ 16( Pi r „-l ’ X n-1 ) + 5<f>( C-2-^-2 )] ■ 

for n = 2,3,---,N—1. 

Replacing n by n + 2 , the standard form becomes 

X n + 3 = X n+2 + [ 23 <Pn + 2 ~ 16( Pn + l + 5( Pn \ ( 2 6 ) 

Again, the local truncation error can be shown to be 

T nJ h) = ^x ,4> ( fu„ )h 3 , for n n e (r n _ 2 ,r n+1 ) [7], 

Results from this method shall be displayed along with the results of the exact solutions for our numerical examples. 
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Adams-Bashforth four-step explicit method: 

x 0 = a 0 , x, = a,, x 2 = a 2 , x 3 = a 3 

h 

X n+i = + ^[ 55< P( r n’ x n )~59<p( r n _„x n _j) + 37& r n _ 2 ,x n _ 2 )-9cp{ r n _ 3 ,x n _ 3 )] , 

where n = 3,4,---,N -1. 

Replacing n by « + 3,wehave 

X „ + 4 = X „ + 3 + ^ [ 55( Pn+3 ~ 59 <Pn + 2 + 32 <P n+ l ~ 9( Pn ] ■ (2V) 

The local truncation error is 
251 

h)= — x>5> ( K )h 4 , for some p n e (r n _ 3 ,r n+I ) [7], 

Results from this method shall be shown along with the results of the exact solutions for the numerical examples. 


Adams-Bashforth five-step explicit method 

x 0 = a 0 , Xj = a j , x 2 = a 2 , x 3 = a 3 , x 4 = a 4 

1901<p( r n ,x n )-2774cp{ r n _ lt x n _, ) + 2616<p( r n _ 2 ,x n _ 2 )- 
1271<p( r n _ 3 ,x n _ 3 ) + 251<p( r n _ 4 ,x n „ 4 ) 

where n = 4,5,---,N — 1. 

Replacing n by n + 4,wehave 

x n+5 = x n+4 +^[l901(p n+4 —2774rp n+3 +2616(p n+2 -1271<p n+1 +251cp n ]. 

The local truncation error is 
95 

T n + i( h )=^Z x< ’(Mnfi ’ for some ju n e ( r n _ 4 ,r n+1 ) [7], 


720 


Adams -Moult on method 

To derive Adams-Moulton implicit k-step method, we can form the backward-difference polynomial P k ( r) through 


( r n+ 1 ’<P( r n + i’X(r n+ 1 ))),(r n , , p(r n ,x(r n ))), •••, (r n+I _ k ,<p( r n+ 1 _ k ,x( r n+I _ k ))) .Since P( r ) is an interpolating polynomial of degree 
k , then for some in (r n+1 _ k ,r n+1 ) , we have. 


cp( r,x( r ) = P k ( r) + 


'(L' X (U) 

(k + 1)! 


(r~r n+I Xr-r n )—( r-r n+1 _ k ). 


Introducing the substitution r = r n + sh , and with dr = hds into P k and with the error term, it implies that 


J <p(r,x(r ))dr = J "' ^ ( -1 f 


v ”y r n ,x(r n ))dr 
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t 


)), 


+ 1 ' (k + 'j' ) r < r - r "+' )( r ~ r n )-( r - r n +I - k )dr 


= Y v '"<p( r „’ x ( r „ )M-i r[ 


r' 

M 

Jo 

[mj 


ds 


h k 


(k + 1)! 


7 \ 0 ( s - 1 > s ( s + 1 )■■■( s + k - 1 ') ( P lk+n (%n’ x ( 4 ))ds. 

ds for various values of m can be evaluated easily as displayed below: 


( s 

-l)s 


f \ 

f J 

~ s 

Jo 

K m ) 


For, 


* = JL 




J, 

' =!:(-! T J ( 


1 [ -s 


m 


f' 

0 1 m 


12 

ds = —- 

J 24 

\s = - 19 


720 


= 1 :(-l ) m \\ 5 Us = - — • 
Jo I m f ir,r 


160 


As a consequence of these evaluations, 
| rp( r, x( r ) )dr = h 

l.k+2 


(Pi r n+l ,x( r n+1 ))--V<p( r n+I ,x( r n+1 )]-—V <Pi r n+1 ,x( r n+I )) + ■ 


-\ 1 g (S~ 1 ) S ( S+1 )'"( S + 1 - k ) < P <k+1> ( Zn’ X ( 4 )) ds - 


(k+1)! 

The coefficient (s — 1 )s( s + 1)(s + 2 )■■■( s + k-1) does not change sign on the interval \0,l] [12]. The weighted mean value 
theorem for integrals can be used to deduce that for a number ju n , where r n+I k < ]u n < r n+1 , the error terrain equation above 


becomes 


-\ 0 ( s-i )s( S + 1 )■■■( s+!-k )<p n+ll ( % n ,x( Z„))ds 

V (s _ 1)s(s + 1) ... (s + k _ 1)ds , 

ik + ni J » 


(k+1)! 

U< k +2), 0 (k+l) 

7 k + 1)! 

or the error term can be written as 


[13]. 


r' 

M 

Jo 

UJ 


ds 


Adams-Moulton two-step implicit method: 

Xg = CCg, Xj =a I 

X n + 1 = X „+Y 5 <P( r n + l ’ X n + l) + 8( Pi r n ’ X n)- ( p{r n ,X n )\ 

where n = l,2,---,N — l. 

To put this in standard form, we replace n by n + 1. Thus, 

X n+2 = X n + l+J^[ 5< P( r n + 2> X n + 2 ) + 8( P( r n+ l’ X n + l )~ ( P{ r n > X n )]. 

The local truncation error is 
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TnJh) = -^x ,4) ( ju n )h 3 , fo rsorne Mn e f r n _„r n+1 )Ui- 

Adams-Moulton three-step method: 

x 0 = a 0 , Xj = aj , x 2 = a 2 , 
h 

X n + 1 = x n+ — \ 9 <P( r n+i’ x n + 1 ) + 19<P( r n ,x n )~5<p( r n _j,x n _j ) + rp( r n _ 2 ,x n _ 2 )] , 

for n = 2,3,---,N-1. 

Replacing n by /i + 2,wehave 
h 

X n+3 = X „ + 2 + ~[«W r n + 3’ X n + 3 )+l 9 <P( ^2 ■ X n +2 ) ~ 5 r n+ l ’ X „ + l ) + <P( r „’ X n j] 

The local truncation is 

T n+ i( h ) = -^ x<5> (H,)h 4 , for M n ^(r n . 2 ,r n+I )[ll 

Adams-B as hforth four-step explicit method: 

x 0 =a Qf Xj =a n x 2 =a 2 ,x 3 = a 3 , 

h 

X n+i = x n + ^l 251<p(r n*i’ x n + i ) + 646(p( r„,x n )~264(p{ r n _ 1 ,x n _ I )+106(p( r n _ 2 ,x n _ 2 )-19<p(r n _ 3 ,x n _ 3 )\, where 

n=3,4,---,N—l. 

Replacing n by n + 3, we have 
h 

x n+4 = x n+3 + — \251cp(r n+4 , x n+4 ) + 646(p{r n+3 , x „ +3 ) - 264(p(r n+2 , x n+2 ) + 106cp(r n+1 , x n+I )-19<p(r n ,x „)] 

The local truncation error is 

T " +l( h)= ~7^ X<6 >( )h " ’ for SOim 0 ( r n~3’ r n + l ) 

lot) 

Numerical Examples 

We will now solve the following problems using some of the methods discussed in this work and the results displayed in 
tables along with the results of the corresponding exact solutions. 

Example 1: Given that 

x = x — r 2 , x(0 ) = 1, 

we obtain the values of x for r =0.1,0.2,0.3,0.4 by the Euler method, Runge-Kutta method of order four, mid-point mle, 
Adams-Bashforth two-step explicit method, Adams-Bashforth three-step explicit method and Adams-Bashforth four-step 
explicit method and then compare the result with the exact solution x( r) = r~ + 2r + 2 - e to obtain the error. 

Example 2: We solve the initial value problem 

xx = r, x(0) = l, 0<r<0.4, 

using the Euler method, Runge-Kutta method of order four, mid-point mle, Adams-Bashforth two-step explicit method, 
Adams-Bashforth three-step explicit method and Adams-Bashforth four-step explicit method with step size h =0.1 and then 

compare the results with the exact solution x(r )= ~Jr 2 +1 . 

Euler method for Example 1 

Using equation (18) x n+1 = x n +h<p n , 
n = 0: x : = x 0 + hrp 0 

r o = 0’ x o =1 
<Po=<P(0,l ) = 1 
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Xj =1+ 0.1x1 = 1.1 
n = 1: x 2 = x! + h(Pi 

^=0.1, x I =l.l 
rp, = qf 0.1,1.1) = 1.09 
x 2 = 1+0.1x1.09 = 1.209 
n = 2: x 3 = x 2 + hcp 2 

r 2 = 0.2 , x , = 1.209 

<p 2 = <p( 0.2,1.209) = 1.169 
x 3 = 1.209 + 0.1x1.169 = 1.3259 
n = 3: x 4 =x 3 + hcp 2 

r t = 0.3 , x 3 = 1.3259 

rp 3 = <p( 0.3,1.3259 ) = 1.2359 
x 4 = 1.3259 + 0.1x1.2359 = 1.449949 
The result of the Euler method is displayed in TABLET 


Table. 1: Euler’s rule for Example 1 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.100000 

1.104829 

0.004829 

0.2 

1.209000 

1.218597 

0.009597 

0.3 

1.325900 

1.340141 

0.014241 

0.4 

1.449949 

1.468175 

0.018226 
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Runge-Kutta method of order four for Example 1 

Using equation (16) x n+1 =x n +—(k I +2k 2 +2k, +k 4 ), 

'6 

where k t = h(p( r n ,x n ), 

k , = h(p\ r +—,x + — k, 1 
2 2 " 2 ') 

k, = lim r + — ,x + — k, \ 

k 4 =hrp(r n +h,x n +k 3 ) 
n =0 : r 0 = 0, x 0 = 1 
k, =0.1xq>(p,l ) 

(p{0,l) = 1 

k, =0.1x1 =0.1 

k 2 =0.1x<p(0.05,1.05) 

(p(0.05,1.05) = 0.10475 


k 2 =0.1x0.10475 =0.10475 
k 3 = 0. l x (p(0.05.1.052375) 
(p(0.05,1.052375) = 1.049875 
k 3 =0.1x 1.049875 = 0.1049875 
k 4 =0.1x<p(0.1,1.1049875) 
(p(0.1,1.1049875) = 1.0949875 
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k 4 =0.1x1.0949875 = 0.10949875 x= l + -[0.1 + 2 (0.10475)+ 2(0.1049875)+ 0.10949875] = 1.104829 

n= 1: r, =0.1, x,= 1.104829 

k, = 0.1x^(0.1,1.104829) 
rp(0.1,1.104829) = 1.094829 
k, =0.1x1.104829 = 0.109483 
k 2 = 0.1x^(0.15,1.159571) 
rp(0.15,1.159571) =1.137071 
k 2 = 0.1x1.137071 = 0.113707 
k 3 =0.1xrp(0.15,1.161683) 
rp(0.15,1.161183) = 1.139183 
k 3 =0.1x1.139183=0.113918 
k 4 = 0.1xrp(0.2,1.218747) 

<p(0.2,1.218747) = 1.178747 
k 4 =0.1x1.178747 = 0.117875 

x 2 = 1.104829+ ^[0.109483+ 2(0.113707) +2(0.113918) + 0.117875] = 1.218597 

n =2: r 2 =0.2, x 2 = 1.218597 
kj =0.1x^(0.2,1.218597) 

<p(0.2,l..218597) = 1.178597 
k, =0.1x1.178597 = 0.117860 
k 2 =0.1x(p(0.25,1.277527) 

<p(0.25,1.277527) = 1.215027 
k 2 =0.1x1.215027 = 0.121503 
k 3 =0.1x(p(0.25,1.279349) 
rp(0.25,1.279349) = 1.216849 
k 3 =0.1x1.216849 =0.121685 
k 4 =0.1x<p(0.3,1.340282) 
cp(0.3,1.340282) = 1.250282 
k 4 =0.1x1.250282 = 0.125028 

x 3 = 1.218597 + ^\ 0 - 117860 + 2 (' 0 ■ 121503)+ 2 (0.121685) +0.125028] = 1.340141 

n =3: r,=0.3, x 3 = 1.340141 

kj =0.1x<p(0.3,1.340141) 

<p(0.3,1.340141) = 1.250141 
kj = 0.1 x 1.250141 = 0.125014 
k 2 =0.1x<p(0.35,1.402648) 
rp(0.35,1.402648) = 1.280148 
k 2 =0. lx 1.280148 = 0.128015 
k 3 =0.1xrp(0.35,1.404149) 
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<p(0.35,1.404149) = 1.281649 

k 3 =0.1x1.281649 =0.128165 

k 4 = 0.1x(p(0.4,1.468306) 

rp ( 0.4,1.468306 ) = 1.308306 

k 4 =0.1x1.308306 = 0.130831 

x 4 = 1.340141 +1[< 1.125014 +2(0.128015 ) +2 ( 0.128165)+0.130831 ] = 1.468175 
The result of the Runge-Kutta method is displayed on TABLE2. 


Table.2: Runge-Kutta method of order four for Example 1 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.104829 

1.104829 

0.000000 

0.2 

1.218597 

1.218597 

0.000000 

0.3 

1.340141 

1.340141 

0.000000 

0.4 

1.468175 

1.468175 

0.000000 
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Mid-point rule for Example 1 

Using equation (19) x n+2 = x n +2h(p n+1 , 

n = 0: x 2 = x 0 + 2hcpj 

r 0 =0, x 0 = 1 , rj=0.1, Xj = 1.104829 

cp,=(p (ft 1,1.104829) = 1.094829 
x 2 = 7 + 2 x 0.1 x 1.094829 = 1.218966 
n = 1: x 3 = Xj + 2h<p 2 

r 2 = 0.2 , x 2 = 1.218966 
<p 2 =(p{0.2,1.218966) =1.178966 
x 3 = 1.104829 + 2 x 0. lx 1.178966 = 1.340622 
n = 2: x 4 = x 2 + 2hcp 3 

r 3 = 0.3, x 3 = 1.340622 

rp 3 = <p(0.3,1.340622) = 1.250622 
x 4 = 1.218966 + 2 x 0.1x1.250622 = 1.6909 
The result of the mid-point rule is displayed on TABLE 3. 


Table.3. Mid-point rule for Example 1 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.104829 

1.104829 

0.000000 

0.2 

1.218966 

1.218597 

0.000369 

0.3 

1.340622 

1.340141 

0.000481 

0.4 

1.469090 

1.468175 

0.000915 


Adams-Bashforth two-step explicit method for Example 1 

Using equation (25) x n+2 = x ll+1 + ^[3cp n+1 -cp n ], 
h 

n = 0: x 2 =x I + -[3cp I -(p 0 \ 
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n = 1: 


r 0 =0, x 0 = 1, r 3 =0.1, x,= 1.104829 

% = 9(0J) = 1 

(p,=rp ( 0.1,1.104829 ) = 1.094829 

x 2 = 1.104829 + y[3(7. 094829) -l]= 1.219053 

x 3 — x 2 + — [ 3tp 2 — rpj J 

Tj =0.1, x,=l.104829, r 2 =0.2, x 2 =1.219053 

rp, =(p{0.1,1.104829) = 1.094829 
rp 2 = ( 0.2,1.219053 ) = 1.179053 


x 2 = 1.219053+^-[3(1.179053)-1.094829] = 1.341170 

n = 2: x 4 =x 3+ -[39 3 -9 2 ] 

r 2 = 0.2 , x 2 = 1.219053 , r 3 = 0.3 , x, = 1.341170 
rp 2 = 9 ( 0 . 2 ,1.219053) =1.179053 
rp 3 = ( 0.3,1.341170 ) = 1.251170 

x 4 = 1.341170+^-[3(1.251170)-1.179053] = 1.469893 
The result of the Adams-Bashforth two-step method is displayed in TABLE 4. 


Table.4: Adams-Bashforth two-step method for Example 1 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.104829 

1.104829 

0.000000 

0.2 

1.219053 

1.218597 

0.000456 

0.3 

1.341170 

1.340141 

0.001029 

0.4 

1.469893 

1.468175 

0.001718 


Adams-Bashforth three-step explicit method for Example 1 

Using equation (26) x n+3 = x n+2 + — [23tp n+2 -16<p n+1 + 5rp n ] 


n=0: 


n = 1: 


x 3 = x 2 + — \33cp 2 -16CPJ+ 5rp 0 ] 

r 0 =0, x 0 = 1 , r ; =0.1, x, = 1.104829 , r 2 =0.2 , 
x 2 = 1.218597 


% = 9(0J) = 1 

<p, = <p(0.1.1.104829 ) = 1.094829 
<p 2 = <p(0.2,1.218597) = 1.178597 


x 3 = 1.218597 + ^[23(1.178597)-16 (1.094829)+ 5(1)]= 1.340184 
h 

x 4 =x i + J^ [23^ -16rp 2 + 5(p 2 ] 

r ; =0.1, x, = 1.104829 , r 2 =0.2 , x 2 =1.218597 
r 3 =0.3, x 3 = 1.340141 
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<p, = <p(0.1,1.104829) = 1.094829 
<p 2 = <p(0.2,1.218597) = 1.178597 
rp 3 = <p(0.3,1.340184) = 1.250184 

x 4 = 1.340184 + ^[23(1.250184)-16 (1.178597) +5 (1.094829)] = 1.468274 
The result of the Adams-Bashforth three-step explicit method is displayed in TABLE 5. 
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Table.5: Adams-Bashforth three-step explicit method for Example 1 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.104829 

1.104829 

0.000000 

0.2 

1.218597 

1.218597 

0.000000 

0.3 

1.340184 

1.340141 

0.000043 

0.4 

1.468274 

1.468175 

0.000099 


Adams-Bashforth four-step explicit method for Example 1 

h 

Using equation (27) x n+4 = x n+3 + — [55cp n+ 3 -59(p n+2 + 37(p n+ 1 -9(p n ] 
h 

n = 0 : x 4 = x 3 + — \55(p 3 - 59cp 2 + 37(p 2 - 9cp 0 ] 

r 0 = 0, x 0 = 1, r 3 =0.1, x 3 = 1.104829 , r 2 =0.2 , 
x 2 = 1.218597, r 3 =0.3, x 3 = 1.340141 , 

% =<p( 0 ,l) = l 

cp, = rp(0.1,1.104829) = 1.094829 
<p 2 = rp(0.2,1.218597) = 1.178597 
<p 3 = <p(0.3,1.340141) = 1.250141 

x 4 = 1.340141+ < ^-j[55 (1.250141)- 59 (1.178597)+ 37 (1.094829)-9(1)] = 1.468179. 


The result of the Adams-Bashforth four-step explicit method is displayed in TABLE 6. 


Table. 6 : Adams-Bashforth four-step explicit method for Example 1 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.104829 

1.104829 

0.000000 

0.2 

1.218597 

1.218597 

0.000000 

0.3 

1.340141 

1.340141 

0.000000 

0.4 

1.468179 

1.468175 

0.000004 


Euler method for Example 2 

Using equation (18) x n+1 = x n +hcp n , 

n = 0 : x 3 = x 0 + hrp 0 

r 0 =0, Xg = 1 

<Po=<P(0,l) = 0 
Xj = 1 + 0 . 1 x 0 = 1 
n = 1 : x 2 = x 3 + hrpj 

r 3 = 0 . 1 , x 3 = 1 
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(p,=(p{ 0 . 1,1 ) = 0.1 

x 2 = 1 + 0 . 1 x 0.1 = 1.01 
n = 2 : x 3 = x 2 + hrp 2 

r 2 = 0.2 , x 2 = 1.01 
<p 2 =<p( 0.2,1.01) = 0.198020 
x 3 = 1.01 + 0.1x0.198020 =1.029802 
n = 3: x 4 = Xj + h<p 3 

r 3 = 0.3, x 3 = 1.029802 

<p 3 = <p(0.3,1.029802) = 0.291318 

x 4 = 1.029802 + 0.1x0.291318 = 1.058934 
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The result of the Euler method is displayed in TABLE 7. 


Table.7: Euler’s rule for Example 2 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.000000 

1.104988 

0.004988 

0.2 

1.010000 

1.019804 

0.009804 

0.3 

1.029802 

1.044031 

0.014229 

0.4 

1.058934 

1.077033 

0.018099 


Runge-Kutta method of order four for Example 2 

Using equation (16) x n+1 = x lt + — (kj +2k 2 +2k ? +k 4 ), 

6 


where 


r o = 0 ’ 


k, = hrp( r n ,x n ), 

k 2 = h(p[r n +^,x n + ^k^ 

k, = hcp\ r +— ,x + — k, ] 

k 4 =h<p(r n +h,x n +k 3 ) 
x 0 = ! 

k,= 0 . 1 x(p{ 0 ,l) 

(p( 0 ,l )=0 
k, = 0 . 1 x 0 = 0 
k 2 =0.1x(p(0.05,l ) 

<p(0.05,l)=0.05 
k 2 = 0.1x0.05 = 0.005 
k 3 =0.1x(p(0.05,1.0025) 
rp (ft 05,1.0025) = ft 049875 
k 3 = 0.1x0.049875 = 0.004988 
k 4 =0.1xrp(0.1,1.004988) 
rp(0.1,1.004988) = 0.099504 

k 4 = 0.1x0.009950 = 0.009950 x, = 1+ ^[0 + 2(0.005) + 2(0.004988)+0.009950)\ =1.004988 
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n =1: r, =0.1, x 3 = 1.004988 
k, = 0.1x^(0.1,1.004988) 
rp(0.1,1.004988) =0.099504 
k, =0.1x0.099504 = ft 009950 
k 2 = 0.1x^(0.15,1.009963) 

(p (0.15,1.009963) = 0.148520 
k 2 =0.1x0.148520 = ft 014852 
k 3 =0.1x<p(0.15,1.0112414) 
rp(0.15,1.012414) =0.148161 
k 3 =0.1x0.148161=0.014816 
k 4 =0.1x(p(0.2,1.019804) 

<p(0.2,1.019804) =0.196116 
k 4 = 0.1x0.196116 = 0.019612 

x 2 = 1.04988 + ^ [ft 009950 + 2(0.014852) + 2(0.014816) + 0.019612 ] = 7. ft/ 9804 

n =2: r 2 =0.2, x 2 = 1.019804 
k, = 0.1 xrp(0.2,1.019804) 
cp(0.2,1.019804) =0.196116 
k, = 0.1x0.196116 = 0.019612 
k 2 =0.1x<p(0.25,1.029610) 
cp(0.25,1.029610) =0.242810 
k 2 = 0.1x0.242810 = 0.024281 
k 3 =0.1x(p(0.25,1.031945) 
cp(0.25,1.031945) = 0.242261 
k 3 = 0.1x0.242261 =0.024226 
k 4 =0.1x(p(0.3,1.044030) 
cp(0.3,1.044030) =0.287348 
k 4 =0.1x0.287348 = ft 028735 

x 3 = 1.019804 + ^[0.019612 + 2(0.024281) + 2 (0.024226)+0.028735\ = 1.044031 

n =3 r 3 =0.3, x 3 = 1.044031 

k, =0.1x(p(0.3,1.044031) 

<p(0.3,1.044031) =0.287348 
k, =0.1x0.287348 = ft 028735 
k 2 =0.1x(p(0.35,1.058399) 
cp(0.35,1.058399) =0.330688 
k 2 =0.1x0.330688 = ft 033069 
k 3 =0.1x(p(0.35,1.060566) 

<p(0.35,1.060566) = 0.330012 
k 3 = 0.1x0.330012 = 0.033001 
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k 4 =0.1x<p(0.4,1.077032) 
tp(0.4,1.077032) = 0.371391 
k 4 = 0.1x0.371391=0.037139 

x 4 = 1.044031 + - [ 0.028735 + 2 ( 0.033069 ) + 2 ( 0.033001 ) + 0.037139 1 = 1.077033. 
6 L J 

The result of the Runke-Kutta method is displayed in TABLE 8. 


Table.8: Runge-Kutta method of order four for Example 2 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.004988 

1.004988 

0.000000 

0.2 

1.019804 

1.019804 

0.000000 

0.3 

1.044031 

1.044031 

0.000000 

0.4 

1.077033 

1.077033 

0.000000 


Mid-point rule for Example 2 

Using equation (19) x n+2 = x n +2hcp n+I , 

n = 0 : x 2 = x 0 + 2hcp I 

r 0 =0, x 0 = 1 , r t =0.1, Xj = 1.004988 

(p,=<p (ft 1,1.004988) = ft 99504 
x 2 = 1 + 2 x 0.1x0.099504 = 1.019901 
n = 1: x 3 = Xj + 2hcp 2 

r 2 = 0.2 , x 2 = 1.019901 
<p 2 = <p(0.2,1.019901) = 0.196097 
x 3 = 1.004988 + 2 x 0.1x0.196097 =1.044207 
n = 2: x 4 = x 2 + 2hcp 3 

r 3 =0.3, x 3 = 1.044207 

<p 3 = <p(0.3,1.044207) =0.287299 
x 4 = 1.019901 + 2 x 0.1x0.287299 = 1.077361 
The result of the mid-point rule is displayed in TABLE 9. 


Table.9: Mid-point rule for Example 2 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.004988 

1.004988 

0.000000 

0.2 

1.019901 

1.019804 

0.000097 

0.3 

1.044207 

1.044031 

0.000176 

0.4 

1.077361 

1.077033 

0.000328 


Adams-Bashforth two-step explicit method for Example 2 

Using equation (25) x n+2 = x n+I + ^ [3<p n+I -<p n \, 

n = 0: x 2 =x I + ^[3<p 1 -<p 0 ] 

r 0 =0, x 0 —l , r ; =0.1, Xj = 1.004988 
(p 0 =(p[p,i) = 0 
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<p,=<p (ft 1,1.004988 ) = ft 099504 

x 2 = 1.004988 + l ^-[3(0.099504)-0] = 1.019914 

n=l: Xj = x 2 + ^ \3(p 2 - (p, ] 

r 3 =0.1, x,= 1.004988, r 2 =0.2, x 2 = 1.019914 

<p j = <p (ft 1,1.004988) = ft 099504 
rp 2 = (0.2,1.019914) =0.196095 

x 2 = 1.019914 + ^-[3 (0.196095)-0.099504] = 1.044353 

n = 2: x 4 =x 3+ -[3cp 3 -cp 2 ] 

r 2 = 0.2 , x 2 = 1.019914 , r 3 = 0.3 , = 1.044353 

rp 2 = (p (0.2,1.019914) = 0.196095 
<p 3 = (0.3,1.044353) = 0.287259 

x 4 = 1.044353 + [3 (0.287259)-0.196095] = 1.077637 

The result of the Adams-Bashforth two-step explicit method is displayed in Table 10. 


Table.10: Adams-Bashforth two-step explicit method for Example 2 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.004988 

1.004988 

0.000000 

0.2 

1.019914 

1.019804 

0.000110 

0.3 

1.044353 

1.044031 

0.000322 

0.4 

1.077637 

1.077033 

0.000604 


Adams-Bashforth three-step explicit method for Example 2 

Using equation (26) x n+3 = x n+2 + [23cp n+1 -16cp n+1 + 5<p„ ] 

h 

n = 0: x 3 =x 2 + — \23cp 2 -16cp I +5(p 0 \ 

r 0 =0 , x 0 = 1 , r 3 =0.1 , x 3 = 1.004988 , r 2 =0.2 , 
x 2 =1.019804 
<Po=<p(0,l) = 0 

(p, =cp(0.1,1.004988) = 0.099504 
rp 2 =rp(0.2,1.01904) = 0.196116 

x 3 = 1.019804 + [23 (ft 196116)-16 (0.099504) + 5 (0)] = 1.044126 

h 

n = l: x 4 =x 3 + — \23cp 3 -16(p 2 +5(p I ) 

r t =0.1, x 3 =1.004988 , r 2 =0.2 , x 2 = 1.019804 

r 3 =0.3, x 3 = 1.044126 

(p, = rp (ft 1,1.004988) = ft 099504 

<p, = q>(0.2,l.019804) =0.196116 
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rp 3 = <p(0.3,1.044126 ) = 0.287322 

x 4 = 1.044126 + ^ [23 ( 0.287322 ) -16(0.196116) + 5(0.099504 )] = 1.077193 
The result of the Adams-Bashforth three-step method is displayed in TABLE 11. 


Table.11: Adams-Basliforth three-step explicit method for Example 2 


r 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.004988 

1.004988 

0.000000 

0.2 

1.019804 

1.019804 

0.000000 

0.3 

1.044126 

1.044031 

0.000095 

0.4 

1.077193 

1.077033 

0.000160 
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Adams-Bashforth four-step explicit method for Example 2 

h 

Using equation (27) x n+4 = x n+3 +— [55(p n+3 -59<p n+2 + 37(p n+I -9(p n ] 
h 

n = 0: x 4 = x, + —[-5-5<A) ~~ 59cp 2 -\-37cpj- 9cp 0 ] 

r g =0 , x 0 =7, r, =0.1, x, =1.004988 , r 2 =0.2 , 
x 2 = 1.019804 , r 3 =0.3 , x 3 = 1.044031 , 

<Po=(p{0, i) = ° 

(p,=(p(0.1,1.004988) = 0.099504 
<p 2 = <p(0.2,l.019804) =0.196116 
<p 3 = <p(0.3,1.044031) = 0.287348 

x 4 = 1.044031 + ^[55(0.287348)-59(0.196116)+37(0.099504)-9(0y\ = 1.077010 
The result of the Adams-Bashforth four-step method is displayed in TABLE 12. 


Table.12: Adams-Bashforth four-step explicit method for Example 2 


R 

x(r) 

Exact solution 

Error 

0.0 

1.000000 

1.000000 

0.000000 

0.1 

1.004988 

1.004988 

0.000000 

0.2 

1.019804 

1.019804 

0.000000 

0.3 

1.044031 

1.044031 

0.000000 

0.4 

1.077010 

1.077033 

0.000023 


III. CONCLUSION 

In this study, it is seen that the multistep methods are derived using Taylor series expansion and numerical integration. The 
numerical integration approach uses the interpolatory polynomial which is determined by some data points to approximate 
the solution to the differential equations.Here, we employed different single-step and multistep schemes in solving non-stiff 
initial value problems of ordinary differential equations. We found out that,unlike the single-step methods, the multistep 
methods attempt to gain efficiency by using information from all previously computed steps to compute the next solution 
value. From the results of our numerical examples, it is established that the multistep methods, though involving more 
computational effort, clearly show superiority in terms of accuracy compared to the single-step methods. The linear 
multistep methods being discussedalso show stability, and hence ensure convergence.Consistency is seen to hold for the 
multistep methods since they have orders that are greater than or equal to 1 (i.e. A>1). 
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